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Abstract 

This paper presents an alternative proof of the connection between the Ising model on a 
planar graphs G and the set of non-backtracking walks on G first reported in [She60]. The 
techniques used give concise formulas for the spin-spin correlation functions of the Ising model 
without boundary conditions and the Ising model with plus boundary conditions, and as a 
corollary of these formulas the spinor holomorphic fermion observable used in [CHI12] is nat- 
urally identified. The paper is combinatorial and makes use of Viennot's theory of heaps of 
pieces [Vie86]. 

1 Introduction 

The connection between the Ising model and non-backtracking walks began with a paper by Sher- 
man [She60] that made rigorous the Kac-Ward [KW52] approach to the Ising model on 1? . Roughly 
contemporaneous to Sherman's work was that of Vdovichenko [Vdo65, MuslO]. Several papers 
discussing the work of Sherman and Vdovichenko have been published in the intervening years; 
of particular note is an exposition by Bourgoyne [Bur63] . More recently Loebl [Loe04] and Cima- 
soni [CimlO] have obtained generalizations of the expressions of Sherman and Kac-Ward that apply 
to arbitrary finite graphs by embedding them in closed orientable surfaces of sufficiently high genus. 

The new aspects of this paper are that expressions for correlation functions are derived in 
addition to expressions for the free energy, that the paper treats all planar graphs on equal footing, 
and that the proof utilizes established tools of combinatorics, namely the theory of heaps of pieces 
introduced by Viennot [Vie86]. Moreover, the expressions derived for the correlations naturally 
identify a key component of recent works on the Ising model, the spinor holomorphic fermion used 
by [CHI12] . The rest of the introduction summarizes the results of the paper, leaving proofs for 
the remaining sections. No prior knowledge of the Ising model is needed to understand the main 
result of the paper. 

A graph whose vertices all have even degree is called even. The central result of this paper is 
a new proof of the fact that the generating function for even subgraphs of a finite planar graph G 
can be expressed as the exponential of a sum of weighted non-backtracking walks in the graph G. 
Some terminology will be needed. A walk 7 = (71,. . . , 7| 7 |+i) is a sequence of adjacent vertices in 
a graph, and I7I is the length of 7. A walk is non-backtracking if jj 7^ 7j+2 f° r 1 < j < M — 2. Let 
r n b(.ff, x, y) denote the set of non-backtracking walks on a planar graph H that begin at x and end 
at y, and T n b(H) = U x r n b(.ff, x, x). For x 7^ z let Z(x,y,z) denote the angle in (— it,tt) between 
the vectors y — x and z — y when H is embedded in the plane; see Figure 1. 
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Figure 1: The turning angle 6 = Z(x, y, z) is the exterior angle of the polygonal segment (x, y, z). 



Define a weight on non-backtracking walks by 

| 7 | 

, r . z (7,',7,'+l,7,- + 2) 

= -l[K ljlj+1 e* J ^ , (1) 

3=1 

where {K xy \ xy £ E(G)} is a collection of indeterminates and i is the imaginary unit. In the 
definition of ^(7) the indices of the vertices jj are taken mod I7I, i.e., 7| 7 |+j = 7j- 

Theorem. The generating function of even subgraphs of a finite planar graph G can be expressed 
as 

E n pf E A P) 

HcG xy£E(H) \76r nb (G) 
_ffeven 

The Ising model on a graph G = (V,E) is the probability measure (-)l on configurations 
a G {±l}l y l given by 

F (°") = ~z II ex P( L xyO- x o-y)/Z, Z = ^ II exp(L xy a x o-y). (3) 

xy£E ae {±i}\V\ xy£E 

The value a x for x £ V(G) is called the spin at the vertex x, and Z is called the partition function 
for the Ising model on the graph G. The quantities L xy £ C are called couplings. There is another 
closely related model, called the Ising model with plus boundary conditions. This is the measure 
(') L,pius defined by the same equations, but with the sums restricted to a x = 1 for a specified set 
dV C V of boundary vertices. The connection between the Ising model on a graph G and even 
subgraphs of G is established by the so-called high temperature expansion: 

Proposition 1.1. The partition function Z of the Ising model on the graph G can be expressed as 

Z = 2 |y| Yl coshi^e ^2 Yl tann ^y> ( 4 ) 

eS-B HcG xyeE(H) 

H even 

where H even means deg(v) is even for all v € V(H). 

Proof. See [Bax82]. □ 

Hence, for applications to the Ising model the indeterminates K xy will be set to be equal to 
t&nhL xy . An immediate corollary of Theorem 1 and Proposition 1.1 is an expression for the log of 
the partition function as a sum over weighted non-backtracking walks; this quantity is called the 



free energy of the Ising model, and is the quantity that previous combinatorial approaches to the 
Ising model have discussed. Two further corollaries of Theorem 1 are expressions for the spin-spin 
correlation function for both the Ising model and the Ising model with plus boundary conditions 
on planar graphs; for the rest of this section assume G is planar. The precise statements of the 
corollaries and the definition of the winding number of a non-closed walk are given in Section 5. 
Let W( 7 , z*) denote the index (or winding number) of a path 7 about z* G V(G*), where G* is the 
graph dual to G. For x*, y* G V(G*) a walk 7 G T n b is called x*,y* odd if VF(7,x*) + W(7, y*) is 
odd. 

Corollary 1.2. Let G be a finite planar graph and x*, y* G V{G*) denote dual vertices corresponding 
to an arbitrary choice of faces adjacent to x,y G V(G). The spin-spin correlation (a x a y )L for the 
Ising model on G with couplings L = {L xy } is 

(a x a y ) L = (a^/V>l E ^(7)(-1) W(t ^ )+W(7 ' !/ * ) , (5) 

■y&r nh (G,x,y) 

where K xy = t&nhL xy , w^j) = c a j ) w('y), where c a b is an explicitly computable unit modulus constant, 
and 

(At**/V>£ = exp(- 2^ tt)- (6) 

7Gr nb (G) 171 

x*,y* odd 

Remark 1.3. The quantity {nx*l l -y*)L is defined by Equation (6), and the symbols fi x have no 
independent definition. This notation has been chosen to make a link with the physics literature 
on the Ising model. No knowledge of this connection is needed to understand this paper. 

Let G[ ow denote the graph dual to G, with boundary dV(G\ ow ) the dual vertex v* corresponding 
to the outer face of G. 

Corollary 1.4. Let G be a finite planar graph, and consider the Ising model with plus boundary 
conditions on G\ ow . If the couplings L x *y* are Kramers- Wannier dual to the couplings L xy on G, 
then the spin-spin correlation for spins at vertices x*,y* is given by 

1 \ , ^(7), ._. 

K*cv) L * )P i us = exp(- 2^ ~up' 

x* ,y* odd 

where K xy = tanh.L xy and w^) = c a i,w{^), with c a b an explicitly determinable constant of unit 
modulus. 



Combining the corollaries yields 
{o- x a y )L 



£ ^( 7 )(-1) w (^*)+ w (t^). (8) 



Recent work of Chelkak, Hongler, and Izyurov has shown the conformal covariance of multispin 
correlations in the Ising model at criticality, using what they term spinor holomorphic fermion 
observables [CI11, CHI12]. Equation (8) identifies the observable in the case of two spins with the 
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weighted sum of non-backtracking walks on the right-hand side of the equation, which is a different 
representation of the observable than appears in the literature. 

During the preparation of this work the author learned that Corollaries 1.2 and 1.4 on correlation 
functions have been found independently by Kager, Lis, and Meester [KLM12]. In addition, these 
authors have obtained convergence results for the correlation expansions as part of a program to 
provide combinatorial proofs of results about the Ising model that currently rely on the Onsager 
solution [Ons44]. 

1.1 Preliminaries and Outline of Paper 

A graph will mean a finite graph without loops or multiple edges unless otherwise noted. G will 
denote a fixed background graph, and |V(G)| and |£^(G)| will denote the number of vertices and 
edges in G. G* = (V*, E*) will denote the graph dual to G when G is planar and embedded. Edges 
{x, y} will be abbreviated xy when it will not cause confusion. It will be assumed that each vertex 
v of the graph carries a cyclic order on the edges incident to v, and that there is an associated angle 
Z (ej,ej+{) = — Z (ej+i, ej) <G {—it, 0] such that the sum of the angles at a fixed vertex is 2tt, i.e., 

d v 

£z( ej ,e j+ i) = 27r, (9) 

3=1 

for all v € V(G), where e^+i = e\ and d v denotes the degree of the vertex v. The angle Z (ej,ej + k) 
is defined to be Yli=j _1 ^ ( e ^> e i+i)- 

Remark 1.5. The angles defined above can be thought of as a type of local embedding of the 
graph into M 2 at each vertex. This enables turning angles to be assigned to matchings, which are 
introduced in Section 2.1. 

The outline of the rest of the paper is as follows: 

• Section 2 introduces loops and shows that the generating function for even subgraphs is the 
same as a generating function for edge disjoint collections of loops. 

• In Section 3 geometric properties of planar curves are used to factorize the weight on collec- 
tions of loops, and it is shown that the factorized form of the weight is local. 

• Section 4 uses the theory of heaps of pieces, along with a combinatorial loop-erasure argument, 
to show that the generating function of collections of edge disjoint loops is equivalent to the 
generating function of non-backtracking walks beginning and ending at the same vertex. 

• In Section 5 the proof of the main theorem is completed, and corollaries for the Ising model 
are presented. 

Remark 1.6. The Ising model as defined in this paper is occasionally called the Ising model with 
free boundary conditions in the literature. Throughout this paper the Ising model with no qualifiers 
will refer to the measure defined by Equation (3); for grammatical reasons it will also be called the 
Ising model without boundary. 
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2 Decomposing Even Graphs 



The main idea in this section is to represent even subgraphs in terms of orientable objects. This 
is done by choosing, for each vertex v of the subgraph, a matching of the edges incident to v. A 
choice of matchings decomposes an even subgraph into a collection of closed loops, each loop being 
formed by picking a sequence of edges that are matched to one another. 

Section 2.1 constructs a weight on matchings that such that the sum of the weight of matchings 
corresponding to a given even subgraphs is equal to the weight of the subgraph itself. Section 2.2 
shows that matchings are equivalent to loop decompositions of even subgraphs and then transfers 
the weight on matchings to a weight on collections of loops. 

2.1 Matchings of Even Graphs 

Definition 2.1. The line graph C(H) of a graph H is the graph with vertices E(H) and edges 
{{xy,yz} | xy,yz G E(H)}. 

The vertices of the line graph will be called half-edges. For each vertex v G H let l(v, H) denote 
the subgraph induced in C(H) by the set of vertices {vw G E(H)}. Figure 2 illustrates an example. 
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Figure 2: Part of the line graph C{H) superimposed on H, for H = I?. The subgraph l(0,H), which is 
not properly embedded in K 2 , is comprised of the square vertices and dashed edges. 

Definition 2.2. Two edges e\e2, /1/2 € t(v, H) cross at v if e\ < f\ < e2 < /2 in the cyclic order on 
edges at v. 

Remark 2.3. The definition of a crossing has a geometric motivation. If H is embedded in the 
plane, i(v , H) is naturally (possibly non-properly) embedded in the plane by placing the vertices 
of l(v, H) at the midpoints of the edges of H. If the cyclic orientation on edges arises from the 
orientation of the plane, two edges cross if and only if they intersect transversely as line segments 
in the plane. 

Definition 2.4. Let H = (V,E). A perfect matching (or dimer cover) of H is a collection of edges 
m C E such that for any v G V there exists a unique e G m with v G e. 

Define the weight w c {m) of a matching on i(t> , by 

w c (m) = (-l) cr M v^^;, (10) 

where cr(m) denotes the number of edges in the matching that cross at v. 
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Lemma 2.5. Let H be an even subgraph of G, and v £ V(H). Then 

w c {m) = 

II ^ ( n ) 

meM(t{v,H)) xyeV(t(v,H)) 

Proof. Choose an edge incident to the vertex v to get a linear order on vertices adjacent to v from 
the cyclic order on edges adjacent to v. Let A XV)VV = ±y/K xv K yv for xv,yv in E(H), with + 
chosen if x > y in the linear order and — otherwise. Set A xv>xv = 0. With these definitions the 
left-hand side of Equation (11) is the Pfaffian of the skew symmetric matrix A. Using the fact that 
the Pfaffian of (Cj^kBjk) equals f\ €j times the Pfaffian of the matrix B and that the Pfaffian of the 
skew-symmetric matrix that is all 1 above the diagonal is 1 proves the claim. See [Ste90] for proofs 
these facts as well as the fact that the left-hand side of Equation (11) is the Pfaffian of A. □ 

Extend the weight w c on matchings of subgraphs l(v, H) to a weight on sets of matchings 
multiplicatively, i.e., w c {{m Vj }j =1 ) = IJjLi w c {m Vj ). Under this weight, the collection of matchings 
corresponding to an even subgraph H have the same weight as the subgraph itself. Formally, 

Proposition 2.6. Let H be an even graph. Then 

H K xy = Yl MM) (12) 

Proof. Let V = V(H). By distributivity 

Y w c ((m v )) = Y 1.1 w c(m v ) (13) 

("»«)en„ e v M(i(v,H)) K)en„ ev M(l(v,h)) vev 

= 11 E Mmv) (14) 

veV m v eM(i(v,H)) 

= n n (15) 

veV xyeV(i(v,H)) 

where the final equality follows from Lemma 2.5. Observe that each half edge xy belongs to only 
the induced subgraphs l(x,H) and i(y,H), so for each half-edge xy a factor \jK xy is contributed 
exactly twice. □ 

2.2 Loops and Perfect Matchings 

A walk 7 in a graph H is a sequence 7 = (71, ... , 7^) of vertices 7^ G V(H) such that 7j7j+i € E(H). 
If 71 = a and 7^ = b then 7 is a walk from a to b in H. The length I7I of a walk is one less than 
the number of vertices in the walk. A walk is non-backtracking if 7j-i 7^ 7j+i for 2 < j < \ j\ and 
is closed if 7| 7 i+i = 71, and is (edge) simple if {7 m ,7m.+i} = {7n>7n+i} implies m = n. 

Let T(H,a,b) denote the collection of all walks from a to b in H. Define T(H,a) = T(H,a,a), 
T(H) = U a T(H,a), and define r n b(iJ, a, b), T n ^(H,a) and so on similarly, with the subscript nb 
indicating the walks are non-backtracking. When the graph H in which the walks take place is 
clear the H in the notation will be omitted. 

Let £l c (H,a) be the set of closed simple walks from a to a in H, and let £l c (H) = U a fl c (H 1 a). 
Define two closed simple walks to be equivalent if they are equivalent as cyclic sequences, or if one 
is the reversal of the other. Denote the set of walks from a to a under this equivalence relation by 
Cl c (a). Elements of O c will be called loops. 
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Figure 3: The left-hand side of the figure indicates a walk 7 6 f2 c . The right hand side is the collection of 
perfect matchings ^(7) associated to 7. In the matchings the dashed edges belong to a copy of K4, while 
the dotted edges are all matchings on separate copies of K 2 ■ 

Definition 2.7. Two subgraphs H\,H2 C G are edge disjoint, written H\ fl e H2 = 0, if E(H\) PI 
E(H2) = 0. Two subgraphs which are not edge disjoint are said to edge intersect. {Hi, . . . , -fffc}n e 
will denote a set of k pairwise edge disjoint subgraphs. 

The notion of edge intersection extends to loops, as loops are just subgraphs with additional 
structure. 

Definition 2.8. Let H be an even graph. A decomposition of H is a set {Cj} of edge disjoint loops 
such that UjE(Cj) = E(H). Let T>(H) denote the set of decompositions of an even graph H. 

When H is an even graph there is a bijection <p: T>(H) — > Ylv<=v(H) -M(i(v,H)). The bijection 
is most usefully described in a somewhat informal manner. Let 7 € ft c be given by 7 = (v±, . . . , v^). 
The walk 7 defines a matching m x at each vertex x occurring in 7 by pairing the adjacent edges 
that pass through x, i.e., pairing VjVj + \ to Vj + \Vj + 2 if Vj + \ = x; see Figure 3. Let 0^(7) denote the 
matching produced. Observe that if r/ is a cyclic shift or reversal of 7 then <j>{eta) = 0(7), so 4> can 
be viewed as a map on the equivalence classes of walks Cl c . 

Define _1 as follows. Pick an edge e = xy in a collection of matchings (mj G rLel^//) -M(i(v, H)). 
The matchings at the vertices of e specify two further edges e\ = xz\,&2 = yz2, and the matchings 
at zi,Z2 specify two further edges, and so forth. Continuing in this manner produces a loop C±. 
If not all of the edges in the matchings are used, repeat this procedure. It is clear that c/> _1 is the 
inverse of <fi. 

Extend <\> to a map on V(H) by 4>{{Ci}) = {(f>(Ci)}. 4>~ x an inverse for this extended map <p as 
well, as it acts on the image of each loop Cj independently. The preceding discussion will be used 
in the future, so its conclusion is recorded in a lemma. 

Lemma 2.9. Let H = (V,E) be an even graph. The collection T>(H) of decompositions of H is in 
bijective correspondence with Y\ v€ y M(l(v, H)). 

Lemma 2.9 induces a weight on decompositions via w({Cj}n e ) = w c {(f){{Ci})); similarly define 
cr({Cj}n e ) = cv(4>({Ci})). Lemma 2.9 shows that to study the generating function of even sub- 
graphs it suffices to study the generating function of decompositions; henceforth this viewpoint will 
be adopted. To begin, observe that resolving a collection of matchings into loops distinguishes two 
types of crossings. 

Definition 2.10. Let Ci , C2 G f^ c , with C\ n e C2 = 0. The self intersection C\ ■ C\ of C\ is given by 

d ■ d = cr(d). (16) 
The mutual intersection C\ ■ C2 of edge-disjoint loops C\ and C2 is 

d-C 2 = cr({Ci, C 2 }) - cr(Ci) - cr(C 2 ). (17) 
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Figure 4: The right-hand side of the figure indicates perfect matchings of a collection of induced subgraphs. 
The dashed edges belong to a copy of K4 while the dotted edges belong to distinct copies of K 2 . The left 
hand side is the element of T> given by §~ x . 

Lemma 2.11. Let {Cj}n e £ T^{H) be a decomposition of an even graph. Then 

j<k j 
Proof. As a crossing involves only two edges 

CT ( m {Cj}) = ^2 ( CT ( m Cj u m c k ) - cv(m Cj ) - cr(mcj) + £cr(mc,) (19) 

j<k j 

Lemma 2.11 allows us to rewrite the weight on decompositions in a more geometrically intuitive 
form: 

Corollary 2.12. The weight on decompositions can be written as 

H{C 1 }n e ) = U(- 1 f J ' C3 U(- l f 3 ' Ck II K *v- (2°) 

j j<k xy€UE(Cj) 

Remark 2.13. As mentioned in Remark 2.3, the notion of a crossing has a geometric interpretation 
when H is an embedded planar graph. In this case, a decomposition {Cj}n e joins the vertices of 
each graph l(v, H) along the edges of H, resulting in a set of closed curves in R 2 . Two edges cross if 
and only if they intersect as line segments. Slightly perturbing the positions of the vertices of i(v, H) 
ensures that each intersection is a double point. Henceforth when calculating the intersection of 
walks by viewing them as curves in M 2 it will be assume such a perturbation has been performed. 

3 Factorization of the Loop Weights 

This section presumes that the base graph G is planar, and that the cyclic orientation on edges 
at each vertex agrees with the orientation given by a straight line embedding of G. In other 
words, the background graph is taken to be a planar map whose edges are straight lines. Note 
that such straight line embeddings exist for any planar graph; see [F48]. Further, assume there 
are two specified vertices a 7^ b £ V(G). Utilizing these assumptions allows for significantly more 
convenient expressions for the weight on decompositions to be derived. 
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Vertices of G will occasionally be called primal vertices to distinguish them from the vertices of 
the dual graph G*. A dual vertex x* is said to be adjacent to a primal vertex x if x* is associated 
to a face whose edges contain x. It will be necessary to enlarge the original base graph G into a 
slightly larger graph G'. Choose, for a and 6, adjacent dual vertices a', b', a' ^ b'; these vertices 
will be called the vertices associated to a, b. The vertices of G' are the vertices of G along with 
a', b' , and E(G') = E(G) U {aa' , bb'}. Let a a b denote a vertex simple walk in G* from b' to a'; by an 
abuse of notation a a b will also denote the reverse of this path when the orientation is contextually 
indicated. 

In this section loops will often be denoted by 7, and walks will be represented as sequences of 
edges. It will be necessary to extend walks 7 on G to walks 7' on G'. If 7 € T(G, a, b), define 7' to 
be the walk (a'a, 7, bb') in G'. 

It will also be necessary to specify a point i£K 2 for each point x' such that x is located in the 
face associated to x' and x is not located on any edge or vertex in G' . The point x will be called 
the winding point associated to x'. See Figure 5. 
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Figure 5: The black dotted path is an extended path 7' from a' to b'. The dashed black path is the chosen 
dual path a a b between a' and b that makes extended walks into closed curves. The winding points a and b 
used to define winding numbers are shown as red circles, and the dashed red line is a path a a b as occurs in 
the proof of Proposition 3.11. 

Note that extended walks define a sequence of matchings of half-edges in i(v, G') for v € 7, and 
hence it makes sense, using the terminology of Section 2.1, to speak of the self-intersection of 7', 
as well as the intersection of 7' with another walk rj when 7' and rj are edge disjoint. 

A theorem of Whitney on the double points of plane curves will be used to exploit the assumption 
that G is planar; to state the theorem requires a definition. 

Definition 3.1. The turning number t^j) of a closed smooth curve 7 is the number of rotations 
made by j'(t)/ \ j(t)\ about the unit circle as t moves from to 1. 

Theorem 3.2 (Whitney [Whi37]). Let 7 be a smooth closed curve in M? with at most double points. 
Then 

T ( 7 ) = 7 • 7 - 1 mod 2 (21) 
Definition 3.3. The turning number of a loop is defined to be the turning number of a representative. 
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That is, if a loop 7 has a representative 7 = (71, ... , 7^) G Q C (G) then 



r( 7 ) = 



1 fe 

— ^47j,7 3 -+i) 



2vr . 



7fc+i = 71- (22) 



To see that the definition of the turning number of a loop is well defined it is enough to note 
that it only depends on the adjacent edges in a loop, which coincide for all representatives. The 
sign of result of the sum inside the absolute values on the right-hand side of (22) will depend on 
the representative 7 chosen, but the magnitude will not. 

Remark 3.4. If 7 is a piecewise smooth curve, then one can smooth 7 to a smooth curve and the 
turning number of the resulting curve is equivalent to the turning number obtained by defining the 
rotation of j'(t)/ \~f(t)\ to be the angle between the left and right derivatives at points where these 
two derivatives differ. 

Theorem 3.2 also holds for loops. This is because the number of intersections in a loop is 
equal to the number of intersections of the loop when treated as a curve in M 2 as mentioned in 
Remark 2.13, and the turning number of the loop is equal, mod 2, to the turning number of the 
curves by Remark 3.4. This shows Theorem 3.2 holds as the theorem only depends on the turning 
number of a curve mod 2. 

Proposition 3.5. Let 7 be a loop in a planar graph G. Then 

k 

(-lfM = nexp(yZ( 7 „7, + i)), 7| 7 | + i = 71- (23) 
i=i 

Let W(r], x) denote the index (or winding number) of a closed curve rj about the point x as 
defined by Equation (46) below. The next proposition will be of use in calculating the mutual and 
self- intersection of loops. 

Proposition 3.6. Assume 7 G Q n ^(H,a,b), a' 7^ b' , rj G F^^H), and that 7 and 77 are edge 
disjoint. Then 

(-l)7'-»? = (_l)W(»?,a')+W(»?,&'). ^24) 

Proof. The equality follows from a standard technique for computing the winding number of a 
closed curve 7 about x, namely counting the number of intersections between 7 and any path v 
from x to infinity, provided 7 and v are in general position. Hence, taking v\ to be a curve in IR 2 
from V to infinity that is in general position with respect to rj that connects with 7' and setting 
to be 7' followed by v\ gives 

(_l)W(?7,a')+>V(?7,&') _ ^_^V"2+Vvi (25) 
= (-1)^' □ 

Remark 3.7. In Proposition 3.6 the points a' and b' can be replaced by a and b by the constancy 
of W(r], x) on components of M 2 \ rj as a function of x. 

Corollary 3.8. Let 71,72 G T^H) be edge disjoint. Then 71 • 72 = mod 2. 
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Proof. The mutual intersection is not changed by viewing the walks as being on a graph with extra 
vertices. Introduce two extra vertices v±,V2 at distinct points on a single edge of 71 and consider 71 
to be a path from 17 to vi- The walk 72 intersects this new walk exactly as many times as it does 
71. Applying Proposition 3.6 and observing W(72,i>i) = W (72,172) as edge disjointness implies v[ 
and v' 2 are in the same component of M 2 \ 72 . □ 

Definition 3.9. The loop weight it; (7) is given by 

|7| 

u;(7) = (-l) T W +i n^7i7 J+ i- (27) 
i=i 

Applying Theorem 3.2 and Corollaries 2.12 and 3.8 to Proposition 2.6 gives a representation 
for the generating function of even subgraphs that is the basis for the rest of the paper: 

Theorem 3.10. Let G be a planar graph. Then the generating function of even subgraphs is given 
by 

E II K *y= E EH-w)- ( 28 ) 

HCG xyeE(H) { 7j }n e <=£>(#) j=l 

even 

For the derivation of correlation formulas it will be necessary to have a formula for the self- 
intersection of an extended walk from a' to b'. Recall a a b denotes the vertex simple dual walk from 
b' to a', and let "f'a a b denote the concatenation of 7' with a a b\ this is a closed curve beginning and 
ending at 0! . Another enlarged graph will be needed. By hypothesis G' is a planar map, and hence 
there is an associated embedding of a a b into G*. Let G a be the graph formed by taking the union 
of G' and a a b, with a vertex at each vertex of G' , at each vertex of a a b, and at each intersection 
between the edges of a a b and the edges of G' when the edges are considered line segments in M 2 . 

Proposition 3.11. Let 7 € fi n b>(a, b), a' ^ b' . Then 

(_iyY'-7' = (_\y(~l' a ab)+WW a ab ,a)+W{~i' a ab fi)+l ^9) 

where the right-hand side is computed in the graph G a . 

Proof. The self-intersections of ^' a a b are either self-intersections of 7' with itself or intersections 
between a a b and 7', as a a b is simple. Hence it suffices to compute ^y'a a b • ^y'a a b and subtract 7' • a a b- 
The first quantity is easily computed using Theorem 3.2: 

j'a a b • ia a b = T(-j'a ab ) + 1 mod 2. (30) 

The intersection between 7' and a a b can be computed by noting that a path a a b from a to 6 can 
be chosen such that 7' intersects a a b if and only if ^'a a b intersects a a b, and such that a a b and a a b 
do not intersect. To see how to construct a a b, note that a a b being vertex simple implies that for 
5 small enough the ^-neighbourhood of a a b is a topological disk; take a a b to be a segment of the 
boundary of the ^-neighbourhood that begins at a near a' and ends at b near b'. For 5 small enough 
it is clear that crossing a a b can occur only if a a b is also crossed. See Figure 5. Hence, applying 
Proposition 3.6, 

1 ■ OL ah = ^'a a b ■ OL a b - a a b ■ Ct a b (31) 
= W( 1 'a ab ,a) + W^ , a ab ,b). (32) 

□ 
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4 Totally Nonbacktracking Walks 



The theory of heaps of pieces is used throughout this section, though the ideas are sketched below 
without reference to the theory. A good introduction can be found in either of [Kra06, Vie86]. The 
basic definitions, terminology, and theorems as used in this paper are recalled in Appendix A. 

This section presents a bijection between a class of closed walks and the set of pyramids of 
loops. The correspondence can be informally described as follows: a closed walk 7 is traced, taking 
note of each time an edge e is used more than once. If e* is the last edge used multiple times in the 
walk, then the segment that begins immediately after the second-to-last visit to e* and ends with 
the final visit to e* defines a closed subwalk C\ of 7. If C\ is edge simple, remove it from the path, 
leaving a shorter path 7'. If C\ is not simple, the path 7 is declared invalid. Repeat the above 
procedure for the path 7', removing the identified closed walk C2 if it is simple, and declaring 7 to 
be invalid otherwise. Continuing in this way eventually reduces a closed walk 7 to an edge simple 
closed walk or results in 7 being declared invalid. 

Provided 7 is not declared invalid, the removed closed subwalks {Cj} inductively form a partially 
ordered set by defining a closed subwalk C\ to be greater than C2 if C\ was removed later than 
C2 and shares an edge with C2 , and taking the transitive closure of the relations generated in this 
way. The remaining closed simple subwalk of 7 can then be declared the maximal element of the 
poset, as whenever a closed simple subwalk is removed, it shares an edge with the remaining path. 
The partially ordered set generated by this process is a heap of pieces, and is the output of the 
bijection when given a valid closed walk 7. 

To invert the map, first trace the maximal element. This gives an order on the closed simple 
subwalks that are immediately beneath the maximal element in the partial order. Moreover, it 
specifies an orientation on those elements. Inductively inserting these closed simple subwalks into 
the maximal closed simple subwalk in the given order recovers the original path 7. The rest of this 
section makes the preceding description precise. 

The set of heaps of pieces with piece types A and concurrency relation 1Z will be denoted 
T-L(A,TZ), the trivial heaps by T(A,1Z), and the pyramids by V(A,1Z). Labels for heaps of pieces 
will be denoted by £. ~H'((i c , n e ) and V'(Cl c , n e ) will indicate heaps and pyramids of loops in which 
the maximal elements are given labels in Q c . Walks will be expressed as sequences of edges. By 
convention j k+1 =71 if 7 = (71 , . . . , 7&) € 

Remark 4.1. Heaps of pieces in which the maximal elements are given labels in a different set than 
the rest of the elements are not the objects of study in the theory of heaps of pieces. However, in 
this section they will be a useful structure to encode non-backtracking walks, and in Section 5 they 
will naturally arise after applying the results of heap theory to genuine heaps of pieces. 

Definition 4.2. The loop index L^: T — >■ N is defined by 

L fe (7) = max {k | 3! m > k 3 j(k) = 7(m)} , (33) 

k 

with max0 = —00. 

Define a second map L e by L e {^/) = m, where m is the index specified in the definition of £5(7). 
Define L e {^) = —00 if ^(7) = —00. As the notation suggests, and L e indicate the beginning 
and end of a closed subwalk. 
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Let 7 G T nb , G ?T (fi c , n e ), and define (7, H)\ [a, b) = (7', H') when ( 7a , . . . , 7^) G ft c by 

7' = (7i>---.7o-i,T*,---,7|7|) ( 34 ) 
H' = iJ + ( 7a ,...,76-i)- (35) 

In the definition of if' the walk (7 a , . . . , 7&_i) is considered an element of %*(il c ,n e ) with a 
single piece. The addition defined for heaps of loops with oriented maximal elements is the standard 
addition of heaps of loops, with elements that are maximal in the sum retaining their orientations. 

Definition 4.3. The erasure map e: T n b x "H*($7 c ,n e ) is defined by 

e(7, J ff) = (7,ff)\[L b ( 7 ), J L e (7)) (36) 
Observe that the erasure map fails to be defined on a pair (7, H) if: 

• the loop index map returns a finite value, but (7l 6 ( 7 ), . . . , 7l b ( 7 )-i) ^ so (7, H)\ [L5, L e ) 
is not defined. 

• no edge in 7 is repeated, i.e., the loop index map returns 00. 

The first failure of definition is because the walk that is identified contains a backtracking segment 
- the edge used twice is used in opposite directions. The second failure is because all that remains 
is an edge simple closed walk. It will be convenient to extend the definition of e so that e(7, H) = 
(0, H + 7) when 7 is an edge simple closed walk. 

Remark 4.4. It is not hard to see that there exist 7 G T n b such that 6(7, H) is not defined; see 
Figure 6. 

The fact that the loop erasure map is not always defined naturally singles out the subset of T n b 
on which the erasure map can be repeatedly applied. Let e k (j,H) denote the fe-fold application of 
e to (7,fT), i.e., 

e fe (7)=^^e(7,^), (37) 

k terms 

where by convention e°(7, H) = (7, H). For future convenience define 6^(7) = e fc (7, 0). 

Definition 4.5. A path 7 G r n b is totally non-backtracking, denoted 7 G r tn b, if e fc (7, 0) is defined 
for all k G N such that e fc_1 (7, 0) 7^ (0, H). The index #(7) of 7 G r tnb is the natural number k 
such that e fc ( 7 ,0) = (0,#). 

The loop erasure procedure is now used to define a map from r tn b to V'(Cl c , n e ). Before the 
definition of the map can be given, a few more definitions are needed. 

Definition 4.6. The walk order on E(H) induced by a walk 7 G r n b(-ff) is given by e > 7 / if 
max{A; | 7^ = e} > max{/c | 7fc = /}• 

Recall that a maximal element x in a heap in ?i'(£l c ,Pi e ) has a label in Q c , i.e., £(x) = 
(£(x)i, . . . ,£(x)\£M\) is a closed walk. 

Definition 4.7. Let 7 G r tn b, H G 1~L'(Cl c ,C\ e ). Let x,y be maximal pieces in h. Define y ^ 7 x if 

l{y)x > 1 i{x) 1 . 

As the erasure map removes loops that occur at the end of 7 first, the relation makes loops 
removed later greater than loops removed earlier. 
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Figure 6: The walk 7 illustrated in the figure is a member of T n b, but not of r tn b- The loop erasure map 
cannot remove the loop present as it is not simple. 



Lemma 4.8. Let 7 G r tn b- Then < e k^ totally orders the maximal pieces o/e fc (7, 0) for k < #(7). 

Proof. No two maximal elements have labels that share an edge by the definition of a heap of loops, 
and the label of each maximal element shares an edge with the remaining walk. Hence, as the walk 
order is a total order on the edges that occur in 7, -^ e fc( 7 ) is a total order on the maximal pieces. □ 

Lemma 4.9. Let 7 G T tn b and k < #(7). Let pj denote the piece whose label is the loop identified 
by e 3 . Then p h -< € k^ pj 2 only if j\ < j 2 . 

Proof. The proof is by induction on k. For k = 1 the statement is trivial. Suppose the statement 
holds for k = m. Observe that the removal of the (m + l) st loop does not change the relative order 
of Pi, ■ ■ ■ , Pm, so it suffices to prove p m i -< e fc( 7 ) p m +i, where m' is the maximal j < m + 1 such that 
Pj is a maximal element of the heap in e m+1 (7). 

Note that ^(p m +i)i <e m (y) £(Pm)i, as otherwise the loop erasure map would remove £{p m +i) 
prior to £(p m ). If m' = m + 1 the claim follows, as the initial edges of both £(p m ) and £{p m +i) are 
retained in e m+1 (~f) in the same relative position as in e m (~f), since these two loops cannot share 
any edges. 

The order < gm^J IS cl total order on e m (j), so by transitivity if m! < m then £(p m +i)i < t m (^) 
£(p m ')i- Noting that £(p m+ i) and £(p m ') share no edges, it follows that £(p m +i)i an d £(p m ')i arc 
both contained in e m+1 (j) in the same relative position as in e m (7), so the claim holds. □ 

Given 7 G T tn b Lemma 4.9 defines max if, the maximal piece under ~< e fc( 7 ) of the maximal 
pieces of the heap in e fc (7, 0) for any k < #(7). 

Proposition 4.10. //e # W( 7 ) = (0,F) ; then H is a pyramid. 

Proof. Let (7', H') = e* ( - 7 ^" 1 (7). Note that the label of each maximal piece in H' contains an edge 
in 7', and 7' is an element of f2 c . It follows that 7' is the unique maximal element in H, so H is a 
pyramid. □ 

Proposition 4.10 defines the map which will be shown to be a bijection, by turning a totally 
non-backtracking walk into a pyramid with the maximal element a closed walk. Next, the inverse 
map is constructed. Just as loop erasure was applied multiple times to reduce a walk to a pyramid, 
the inverse map is defined by repeatedly applying an operation that will be called loop addition. 
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Definition 4.11. Let 7 € r tn b, H € %*(f2 c ,n e ), and assume maxH (under is unique and such 
that oj = £(maxii) n e 7 7^ 0. Let 7^ be the maximal edge in 7 which is contained in oj. Let 

7 + w = (7i,...,7 fe _i,a;i,...,W| C |,7 fe ,...,7| 7 |), (38) 

and let — max be the heap of loops formed by removing the maximal piece max H, and orienting 
the labels of the new maximal elements according to the edge they share with 7 + oj that occurs 
last in 7 + oj. The loop addition map e _1 is 

e _1 (0,if) = (maxiT,#-max#), (39) 
e _1 (7,if) = + oj,H -m&xH), 7 / 0. (40) 

Lemma 4.12. Let-y£T tnh . Then e _1 o e fe+1 (7, 0) = e fc (7, 0) /or all k < #(7). 

Proof. Let e fc+1 (7,-£f) = (j',H'). First assume 7' is not null. By Lemma 4.9 maxii was removed 
by the (k + l)st application of e. Hence e _1 identifies the correct loop. The orientations of the 
maximal pieces of H are the orientations induced by the last edge of each maximal piece that is 
contained in 7, so e _1 correctly orients the maximal elements of H by virtue of identifying max if 
correctly. □ 

Lemma 4.13. Let P £ V(Cl c ,n e ), k < \P\. Then e~ k ($,P) = (<y,H) is well-defined, 7 € r tnb; 
and e o e- k (Q), P) = £- fc+1 (0, P). 

Proof. The claim is immediate for k = 1. Suppose e _fc+1 (0, P) = (7, H). By induction the maximal 
pieces of H have labels that edge intersect 7, so each maximal piece x of H is given a non-zero rank 
by max e {e € £(x)}. Further, these ranks are distinct, as no two maximal pieces can have labels 
with an edge in common. Hence the construction used in applying e _1 is well defined, as \H\ > 
by the hypothesis k < \P\. 

To show that 7 G r tn b it suffices, by induction, to show e o er k = e~ fe+1 . This follows from the 
construction of e _1 and e as £(maxH) is the last closed subwalk in e _fc (0, P). □ 

Definition 4.14. Let 7 £ r tn t,. Define the /oop erasure bijection map e: r tn b — >■ P*(n c , n e ) by 

g( 7 ) = e #(7) (7) . (41) 

By Proposition 4.10 the range of e is the set of pyramids of loops whose maximal elements are 
closed walk as claimed in the definition. Lemmas 4.12 and 4.13 show that e is a bijection. Further, 
it is evident from the constructions given that the set of edges present in a totally non-backtracking 
walk 7, counting multiplicities, is the same as the set of edges in the labels of the pieces of the 
pyramid corresponding to 7, again counting multiplicities. Hence the main result of this section 
has been obtained. 

Theorem 4.15. The map e is a bijection, and moreover, \{k \ 7^ = e}\ = \{x G H \ e € for 
all e £ 7. 

Totally non-backtracking walks 7 have the property that a cyclic shift of 7 or the reversal of 7 is 
still totally non-backtracking. This property is an important ingredient of the proof of Theorem 5.6, 
and is a corollary of the bijection just presented: 
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Corollary 4.16. Let 7 = (71, ■ ■ ■ ,7| 7 |) G rW,- Then any cyclic shift of 7 or the reversal of '7 is 
also a totally non- backtracking walk. 

Proof By Theorem 4.15 7 can be written as a pyramid P with a rooted maximal element. Note 
that there is a piece x e G P corresponding to each edge e in the walk. Picking a walk 77 that is a 
cyclic shift of 7 or the reversal of 7 is equivalent to choosing a particular edge e v that occurs in 7 
and an orientation for that edge. Let x v be the piece of the pyramid such that e <G £(x). 

As discussed in [Vie86] heaps of pieces are equivalent to acyclic labelled directed graphs, where 
two vertices are joined by an edge if and only if their labels intersect. If an edge is directed from 
vertex a to vertex b the piece b lies above a. A pyramid is thus an acyclic labelled directed graph 
with a unique vertex with outdegree zero. 

To prove the corollary, trace rj, and each time a new piece is entered, direct this piece towards 
the previous one. This results in an acyclic directed graph with a unique maximal piece x„, as 
each other piece is preceded by another. The graph is acyclic as if adding an edge produced a 
cycle, the most recently entered piece must already have an incoming edge; but this contradicts the 
piece being previously unvisited. By Theorem 4.15 the resulting pyramid corresponds to a totally 
non-backtracking walk. To see that the non-backtracking walk is r] follows from the definition of 
loop addition and the construction of the pyramid. □ 



5 Main Theorem and Corollaries for the Ising Model 

5.1 Proof of Main Theorem 

In this section walks are given as sequences of edges, unless otherwise noted. Define a weight 
w : T nb ->■ C by 

|7| 

0(7) = (-1) T{7) II K *u = U e M^A7j,w))K lj . (42) 

xyej j=l 

To prove the main result requires lemmas showing the compatibility of the weight on loops 
with the gluing operation, and hence with the bijection between pyramids of loops and totally 
non-backtracking walks. 

Definition 5.1. Suppose 7 1 G T n b, 7 2 G Cl c edge intersect. The gluing of 7 1 and 7 2 is the walk 7 
defined by 

7 = (7i\ ■ ■ ■ , 7 S \ 7i , ■ ■ ■ , 7| 7 2| , 7i+i. ■ ■ ■ > 7|V|) ( 43 ) 

where 7* is the first edge in 7 1 that is contained in the loop 7 2 and (jf, . . . , 772]) is the representative 
of 7 2 with initial edge 7* . 

Remark 5.2. This is just the loop erasure construction of Section 4. 

In the next two lemmas it is assumed that the gluing operation is well defined. 

Lemma 5.3. Let 7 be the gluing 0/71 and 72. Then 

w(j) = u) (71)10 (72). (44) 

Proof. The weight is a product of factors, one for each pair of adjacent edges in the walk. Gluing 
preserves the adjacency structure of edges. □ 
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Lemma 5.4. Let 7 be the gluing of 71 and 72. TTien 

W(7,o) = W(7i,a) + W(72,a) mod 2. (45) 
Proof. The index of a curve r\ about a is given by 

W(7/,a)=/— (46) 

with 77 considered a curve in the complex plane. Splitting the domain of integration 7 into its 
constituent parts 71 and 72 gives the result. □ 

The next lemma is conceptually and computationally helpful as non-backtracking walks are a 
significantly less constrained set of walks than totally non-backtracking walks are. 

Lemma 5.5. The sum of w(^) over all closed walks that are not totally non-backtracking is zero. 
That is, 

7er nb \r tnb 

Proof. It suffices to exhibit an involution T on T n b \ r tn b such that w o T = —w. Suppose the 
loop erasure map e can be applied k times to 7 but not k + 1. Define T(7) = e~ k o R o e k ("f), 
where R is the map that reverses the first non-simple closed subwalk of a closed walk. That is, if 
{ojL b , ■ ■ ■ ,wi e ) is the first closed subwalk of ui which is not simple then 

R(U>) = (Ui, . . . ,LOL b -l,UL e ,UL e -l, ■ ■ ■ ,U Lb ,U Le+1 , . . . ,U\ U \). (48) 

R(oj) is an involution and hence T is an involution. Let uiL b = xy. Then ojl £ = yx, and so the 
turning angle of the reversed non-simple closed subwalk differs from the original by 2ir. Hence 
w(R(uj)) = —w(uj). The factorization property of the weight w given by Lemma 5.3 yields w o T = 
-w. □ 

The central theorem of the paper can now be proven. Recall that u>(7) = —^(7). 

Theorem 5.6. The logarithm of the generating function of even subgraphs of a finite planar graph 
G is given by 

io g y. (-n'^m = E # («) 

Ter(n c ,n e ) 7er tnb 1 n 

Proof. That the left-hand side of Equation (49) is the log of the generating function of even 
subgraphs is the definition of a trivial heap of pieces, combined with Theorem 3.10. Let Zq = 

ET e r(^(G),n e )(-l) |T| *(^)- Note that for xy G E{G) 



\ogZ G 



= log Z G \ xy , (50) 



where G \ xy is the graph G with the edge xy removed. The fundamental theorem of calculus then 
implies 

f Kxy d ~ 

log Zq = I log Z G dK xy + log Z G \ xy , (51) 

JO UK X y 
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where Z G is a function of the dummy variables K xy . 

The derivative of log Z G can be identified with a class of pyramids in V: 

' \o g Z G = --±— V w{C) Y(-lf\w(T). (52) 



dK xy ' K xy Z G 

xydE(C) 

The outer sum is over loops C containing xy, and T G is the set of trivial heaps whose elements do 
not intersect the loop C. Applying Theorem A.l gives 

-^-logZ G = --L £ w(C)w(H). (53) 

a x y x y h&hc 

The maximal pieces in the heaps in 7i G edge intersect C, so the sum in equation (53) can be 
identified with 1/2 the sum over pyramids whose maximal element is C, rooted at x and oriented 
towards y or vice versa. Integrating (53) over K xy gives a prefactor of the iV^ 1 , where N xy is the 
number of times xy is contained in the pyramid. Applying Theorem 4.15 gives 



/' 

J o 



\ og z G dK xy = -otf 1 ' ( 54 ) 



dK xy 2N xy 
ne{xy,yx} 



where Lemma 5.3 has been used to write the weight of a pyramid as the weight of a walk. 

Note that the prefactor (2N xy )~ 1 shows that the sum could be considered to be over totally non- 
backtracking walks up to cyclic shifts and reversals; note that all of the cyclically shifted/reversed 
totally non-backtracking walks are still totally non-backtracking by Corollary 4.16. As there are 
2 |7| ways to root and orient a walk, this implies 

f Kxy d ~ v - wh) 

Jo dK xy ^ 2 1 7 | 

^Bxy 

where the condition on 7 is that either xy or yx is contained in 7, and —w('j) = w (7) has been 
used. Applying this argument to log Z G \ xy by selecting another edge wz, and repeating, yields 

>og ZG = £ # <«> 

7er tnb 1 ' 1 

Applying Lemma 5.5 allows the index of the sum to be expanded to all non-backtracking walks 
and completes the proof. □ 

5.2 Consequences for the Ising Model 

Theorem 5.6 combined with Proposition 1.1, immediately yields a representation of the free energy 
of the Ising model in terms of non-backtracking walks. 

Theorem 5.7. The extrinsic free energy F G = log Z of the Ising model on a finite planar graph G 
is given by 

F G = \V\ log2 + Y logcoshL,, + \ IT" (57) 

xyeE 7er- b (G) 171 
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Remark 5.8. By approximating Z 2 by a finite torus and applying the above formula the Onsager 
solution to the Ising model can be recovered. See [She60] for details. 

The starting point to obtain a formula the correlation function for the Ising model is the high- 
temperature expansion formula for the correlation. In what follows the notation (-)l will denote 
the Ising measure with couplings L = {L yz }, and (-)l,p1us wm denote the Ising measure with plus 
boundary conditions with couplings L. The graph G under consideration will be implicitly identified 
by the set of couplings. 

Proposition 5.9. Let G be a planar graph, a / b € V(G). Then 

( \ 

(<7a(Tb)L = 



He£(G ab ) xyeE(H) 
\ ab<=E(H) ) 



E II L *y) > ( 58 ) 

He£(G) xyeE(H) 



where G a b is the graph G with an edge ab added if one is not already present, in which case L ab = 1. 

Proof. This follows from applying the high-temperature expansion to the numerator and the de- 
nominator of the expression that defines (a a o~b)L- □ 

Recall that 7' is 7 begun at a' and ending at b', a a b is the chosen vertex simple dual path from 
b' to a', and j'a a b is the closed walk formed by concatenating these two walks. 

Definition 5.10. Let 7 € T D \y(G, a, b). The winding x) of 7 about x is defined to be W{~i'a a b, x). 
Define a closed walk 7 to be a, b odd if W(7, a) + W(7, b) is odd. 

Theorem 5.11. The spin-spin correlation function (o- a o~b)L for a ^ b is given by 

(a a a b ) L = ( t i a , f i b ,) L £ «)( 7 )(-l)W(7,a)+W(7,5) j (59) 
jer nh (G,a,b) 



where 1^(7) = c a bw{^'), c a b is an explicitly computable modulus one constant, and 

Hi 
h\ 



{/j, x ',H y ')L = e^p(- E IT)- ( 60 ) 



7er nb (G) 

x,y odd 

Proof. For the proof assume that ab is not already an edge in G. The proof is similar, but easier, 
in the case ab <G E(G). Let G a b be the graph with the edge ab added, and set L ab = 1. Define 

Z W (G)= £ (-1) |T| ^(T). (61) 

TeT(n c (G),n e ) 

Corollary 2.12, Proposition 3.6, Proposition 3.11, and Remark 3.7 imply that Equation (58) can be 
rewritten 

Kab) = z^Jg) zagY (62) 
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where C'oL a b is the concatenation of the walk C, extended to begin at a! and end at b', with the 
fixed walk a a b from b' to a' . G a b \ C is the graph G a b with the edges in C removed, and 

u)( 7 ) = (-l)W(7,a)+W( 7 ,6) w ( 7 ). (63) 

The constant c a b in the theorem statement is defined by w(C'a a b)/w(C'); this is constant as 
a a b is fixed, and allows w(C) to be replaced by w(C) in Equation (62). 

Note that the first multiplicative factor in Equation (62) is of exactly the form handled in the 
proof of Theorem 5.6, with w replaced by w. This yields, using the same analysis as in the proof 
of Theorem 5.6, along with the fact that the index factors as given by Lemma 5.4, that 

Ml= £ ™ (7 )(-l)WHM7,*)|^ (64) 

To see that the final ratio of partition functions is equal to (/U a / [iy ) l observe that the numerator 
and denominator are each the generating function of even subgraphs with different weights, and 
apply the same analysis as before to these generating functions, again using that the index factors 
as given in Lemma 5.4. As w(^) — ^(7) = 2w{^) when 7 is a, b odd, and is zero otherwise, the 
formula follows. □ 

Remark 5.12. In [KC71] the quantities fi x i, called disorder operators were defined, and the quantity 
(Hx' fJ>y') l was defined to be a disorder-disorder correlation. 

To conclude this section a proof of the relationship between the disorder-disorder and spin-spin 
correlations is sketched. This relationship was first proven in [KC71] for the case of G = 1? . 

Recall that given a planar graph G the graph Gi ow is the dual graph G* with boundary the 
dual vertex associated to the external face. Using the low-temperature expansion shows that the 
free energy of the Ising model with plus boundary conditions on Gi ow is a generation function for 
even subgraphs on G. See [Bax82] for the details of the low temperature expansion. 

The next definition was first introduced in a landmark paper by Kramers and Wannier [KW41]. 

Definition 5.13. If {L xy }, {L* y } are couplings on G and G* respectively, the couplings are called 
Kramers- Wannier dual if exp(— 2L x * y *) = tanhL X2/ . 

Theorem 5.14. Let a,b £ V(G), a^b, a 1 / b' . The disorder-disorder and spin-spin correlations 
are dual, meaning 

(Va'Vb')L = (o-a'O-b'} L*,plns, (65) 

when the couplings L x * y * = L* y are Kramers- Wannier dual. 

Proof. The method of proof for the low temperature expansion implies that 

/ \ 

' ]T w{H)- £ w{H) ■ (66) 

\a,6evcn a,bodd / 



(0"a'0"b')i*,plus 



Rewriting the sums over even subgraphs as sums over decompositions and noting that the sign due 
to being a, b even or odd is replicated by a factor 

(_l)W(C*,a)+W(C,b) ( g7 ) 
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for each loop C in a decomposition gives 

U „ \ £TeT(n°,n e ) ^( T ) , , 

\0 a' <?b')L*, plus — ^ , (OoJ 

which is equal to Equation (60) by using Theorem 5.6. □ 

5.3 Identification of the Spinor Holomorphic Fermion Observables 

This short section notes that Equation (62) naturally identifies one of the holomorphic fermionic 
observable used in [CI11, CHI12]. As previously described, (/V/"f/) is what [KC71] defined to be 
the correlation of disorder operators; they showed that 

(/V M>') L = (^a'^b')L*, plus, (69) 

where the couplings {L*} are on Gi ow , and the couplings L x * y * are the couplings obtained by 
applying the Kramers-Wannier duality transformation to the couplings L xy . Define 

K/VW) = ^( 7 )(-l) H/(7 ' a)+W(7,5) , (70) 

7er nb (a,6) 

so that Theorem 5.11 can be compactly re-expressed as 

(cr a Vb)L = {Va^a'O h lJLb')(Va'^b')L (71) 
= (VafJ-aiVbVbi}{&a'(?bi)L*, plus, (72) 

SO 

(craVa'VbVb') = 1 T • (73) 

{O a ' O'b'jL*, plus 

At criticality L* = L, combining this with Equation 2.4 of [CHI12] shows that (up to multiplica- 
tive constants) ((JaUa'&blJ'b') is exactly the holomorphic fermionic observable considered in [CHI12]. 
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A Heaps of Pieces 

See either of [Kra06, Vie86] for an introduction to the theory of heaps of pieces. The notation used 
here is the same, though some of the terminology is slightly different. 

A concurrency relation is a symmetric and reflexive binary relation. Let TZ be a concurrency 
relation on a set £>. The set B will be called the set of piece types. A heap of pieces (H, ■<,£) is a 
triple with I: H —> B and (H, ^) a poset such that 

1. If x, y G H and £(x)7Z£(y) then either x ■< y or y ■< x. 

2. The relation ^ is the transitive closure of the relations from the previous condition. 
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The map £ is called the labelling of the pieces. Given a collection of piece types B and a 
concurrency relation 7Z define 

1. W(B,1Z) to be the set of all heaps of pieces, 

2. T(B,1Z) to be the set of trivial heaps of pieces, i.e., heaps of pieces for which £(x)%£(y) for 
any x,y G H. 

3. V(B,1Z) to be the set of pyramids, i.e., heaps of pieces which contain a unique maximal 
element. 

There is a natural notion of addition (or composition) of two heaps of pieces. If (Hi, <i,li) G 
H(B,K) then define 

(Hi, + (H 2 , < 2 ,£ 2 ) = (H 3 , 4) (74) 

by 

1. H 3 = H 1 UH 2 

2. The partial order ^3 is the transitive closure of 

• v\ ^3 v 2 if vi -<i v 2 , 

• vi ^3 v 2 if vi < 2 v 2 , 

• vi ^3 v 2 if vi G Pi, f 2 G P 2 , and ii{vi)Hi 2 (v 2 ) . 

Intuitively the heap (H 2 , ^ 2 ,£ 2 ) is placed on top of the heap (Hi, ^i,£i). Note that this addition 
is not commutative. 

If M C B then %m (&■> ^) is the se t of heaps of pieces with piece types B whose maximal 
elements have labels in M. 

A weight function w : B — > C naturally extends to heaps of pieces via 

w(H) = Yl w(£(x)). (75) 
xeH 

The main theorem we need from the theory of heaps of pieces is the following. 
Theorem A.l. 

h { ] ~ e™^(-i) |t| ^) ( } 
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